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We compare the inclusive transverse momentum spectra of single pious above p^ = 3 GeV/c mea- 
sured in proton-proton {p-p) collisions at y/a — 21.7 - 23.8 GeV, with next-to-leading order (NLO) 
perturbative QCD (pQCD) predictions using recent parametrizations of the parton densities and 
parton-to-pion fragmentation functions. Although the dependence on the theoretical scales is large, 
the calculations can reproduce the experimental results both in magnitude and shape. Based on 
QQ . the existing data and on a pQCD y^-rescaling of the measured spectra, we provide a practical 

t"*^ ' parametrization of the baseline p-p pion transverse momentum spectrum to be compared to nucleus- 

(^ ' nucleus collisions data at ^s„„ = 22.4 GeV. 
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The study of hadron production at large transverse moinenta (p^ 3> Aqcd ~ 0.2 GeV) in hadronic interactions is 

r^ I a valuable testing ground of the perturbative regime of Quantum Chromodynamics (pQCD), providing information 

Oh' on both the parton distribution functions (PDFs) in the proton, and the parton-to-hadron fragmentation functions 

f-^ . (FFs) [1|. In the last years, a renovated interest in high-p^ hadron production has been driven mainly by studies of 

qj ' "jet quenching" phenomena in high-energy nucleus-nucleus (A-A) collisions [2] as well as of the proton spin structure 

f^ , in polarized p-p collisions [d, IJl- ^^ A-A collisions, the observed large suppression of high-p^ hadron yields compared 

'~~' ■ to (appropriately scaled) p-p cross sections [a, Q - attributed to parton energy loss due to medium-induced gluon 

p^ ] radiation 0, Q - provides valuable information on the transport properties of hot and dense QCD matter J2| . The 

^ ■ energy density at which such "jet quenching" phenomena sets in in A-A collisions can signal the possible transition 

(Nl from a hadronic to a deconfined quark-gluon system. Whereas unambiguous signals of high-p^ hadron suppression 

VjO have been found at the Relativistic Heavy-Ion Collider (RHIC) in central Au-Au collisions at y^s^ = 200 GeV [3, [y| 

Cn and 62.4 GeV [9|, Il0l|. o ne cannot draw any firm conclusion yet at the Super Proton Synchrotron (SPS) energies 

^~~. i\/^NN = l'''-3 GeV) [III due to the lack of a vaHd (experimental and/or theoretical) proton-proton reference p^ . 

t^^ ' Recently, the PHENIX collaboration has presented results on high-p^ neutral pion production in Cu-Cu collisions at 

^D ] energies, ^s„„ = 22.4 GeV, close to the SPS range [Tj]. We present here an experimental and theoretical study of 

25 ' the pion p^ -spectrum in p-p collisions required in order to determine the associated nuclear "suppression factor", 

Raa{Pt) o^ {dNAA/dpT)/{dNpp/dp^) in A-A collisions at this center-of-mass (cm.) energy. 



I. INTRODUCTION 



In Section roi we compile and examine all existing experimental spectra for 7r° [Tjj Ej Hi UM, Hi Hi Hi Hi 
and TT^ [2l|, H| at cm. energies in the range ^/s = 21.7 - 23.8 GeV. We notice that most of the data appear 
d ' to be consistent with each other within uncertainties, despite some spread. In Section lllli we compare these data 
to pQCD calculations at Next-to-Leading Order (NLO) accuracy, as implemented in the Monte Carlo programme 
INCNLO j23l Hi . We discuss in some detail the improvements in the model predictions thanks to the use of recent 
FFs [25|. For a choice of renormalization- factorization scales in the low side (n/pj. = 1/3 — 1/2), the calculations 
can reproduce the experimental results both in magnitude and shape within the uncertainties associated with the 
limited knowledge of the parton-to-pion fragmentation functions (FFs) and parton distribution functions (PDFs) 
in this kinematic range. Finally, a practical parametrization of the p-p pion transverse momentum spectrum at 
^/s = 22 A GeV is provided in Section llVI for use as denominator in the determination of the corresponding nuclear 
modification factor in A-A collisions in the low range of energies accessible at the RHIC collider. 
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II. INCLUSIVE PION SPECTRA IN p-p COLLISIONS AT y/s ; 

MEASUREMENTS 



22.4 GEV: EXPERIMENTAL 



Table [J compiles the 13 measurements found in the literature for neutral [Tj, [T^, [T^, [13, HM 13, 123, IIH and 
charged [2l|, ^^ pion production at cm. energies around ^/s = 22.4 GeV at mid-rapidity {y = 0, corresponding 
to laboratory angles 0iab ~ 1- rad in fixed-target kinematics). The data were measured in the 70's at the CERN 
Intersecting Storage Rings (ISR) collider as well as in the 80 's in various CERN and Fermilab (FNAL) fixed-target 
experiments. The corresponding data points (adding to a total of ^220) have been obtained from the Durham 
database [26|]. Assuming isospin symmetry, the 7r° yield is the same as the (7r+ + n^)/2 yield, and thus we can 
use both data sets to get a combined pion reference spectrum. The last column of Table |T] collects the propagated 
experimental uncertainties of the measurements as reported in the original publications. Two types of errors are 
often quoted: (i) those related to energy scale {p^) uncertainties, and (ii) additional systematic and/or absolute 
normalization (usually luminosity) errors. The p^-scale uncertainties have been transformed into an associated 
absolute cross section uncertainty assuming a local power-law distribution with exponent n w 10. We have 
conservatively added all quoted uncertainties in quadrature with the point-to-point errors. We note that at variance 
with the tt" spectra measured at ^/s « 63 GeV [23| , there is no need to account for possible direct-7 contaminations 
in the oldest "non-resolved" pion spectra since, at the lower cm. energies considered here, the prompt photon con- 
tributions start to be significant only above the momentum range {p^ > 6 GeV/c) actually reached in the experiments. 
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TABLE I: Compilation of inclusive pion production data in p-p collisions around y^s — 22.4 GeV and midrapidity: collision, 
center-of-mass energy, piat (for fixed-target experiments), collaboration/experiment name, bibliographical reference, measured 
Pj, range, total number of data points, and associated systematic uncertainties in the measured cross sections. 

Figure [1] shows all the measured pion p^ spectra. The full range of cross sections covers more than 12 decades. 
Unlike with what was observed at ^/s « 63 GeV [27|, the data taken by the various experiments appear in general 
quite compatible with each other both in shape and absolute cross sections, within the experimental uncertainties and 
within the differences expected (at high-p^) due to the slightly dissimilar cm. energies of the various measurements 
(see section HV Ap . The spectra are characterized by an exponential distribution (with inverse slope ~150 MeV) at 
low-py (Py < 1 GeV/c), followed by a power-law with exponent ~ 10, and then a drop at the highest Pj,'s when 
running out of phase-space for particle production, approaching the kinematical hmit {p™'^^ = \/s/2 = 11-12 GeV/c 
at midrapidity). 



III. INCLUSIVE PION SPECTRA IN p-p COLLISIONS AT ^s ; 

CALCULATIONS 



22.4 GEV: NLO PQCD 



The inclusive cross section for the production of a single pion, differential in transverse momentum p^ and rapidity 
y, takes the following form at next-to-leading order (NLO) [24] : 
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FIG. 1: Compilation of all pion transverse spectra measured in p-p collisions in the range y^ = 21.7 - 23. i 
for details). 
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Fi/p are the parton distribution functions (PDFs) of the incoming protons p, D'^(z, fijj) are the parton-to-pion 
fragmentation functions (FFs) describing the transition of the parton k into a pion, and d$ij^k/dp^dy is the 
Born cross section of the subprocess i + j -^ k + X, and K^jk is the corresponding higher-order term (the full 
kinematic dependence is omitted for clarity). In this paper, we use the INCNLO program [23| to compute the 
cross sections, supplemented with various PDFs and FFs sets (see below). The truncation of the perturbative 



series at next-to- leading order accuracy in as, introduces an artificial dependence, with magnitude O i 



of the 



cross section on initial- and final-state factorization scales, fj,F and fijj, as well as on the renormalization scale 
IdR. The choice of scales is to a large extent arbitrary. One often uses as "standard" choice the hard scale of 
the process, e.g. /i/j = ijf = M ff = Pt- ^ more theoretically sound solution is given by using the Principle of 
Minimum Sensitivity (PMS) [2^. Phenomenological comparisons of pQCD results at various orders (LO, NLO, 
NNLO) among each other and against various experimental data sets (for charm and beauty, top, Z, W bosons, ...) 
indicate that choosing a relatively low range of scales n/p^ = 1/3-1/2 provides effectively a reduced sensitivity to 
higher-order effects [29|. We thus use fin — fip — ^g = Pr^/n, with variation between n — 2-3. At small p^ and 
for the scale p^/i, the factorization scale approaches the starting scale Qo of the PDF evolution, where the parton 
densities are not constrained by data. To avoid this problem, we only compute the pion spectra above^ p^ = 3 GeV/c. 



^ Whenever it becomes smaller than the minimum Qq allowed by the PDF or FF parametrization, the hard scale Q is frozen at Qq. 
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FIG. 2: Comparison of pion transverse spectra measured in p-p collisions at -y/s ^ 21.7-23.8 GeV to NLO pQCD predictions. 
The left (right) plots are for theoretical scales /i = p^/3 {p-r/2). Two sets of PDFs (MRST04 and GTEQ6.1M) and three 
FFs (AKK05, AKK08, DSS, from top to bottom) are used. 
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FIG. 3: Ratio of data over pQCD for pion transverse spectra in p-p collisions at ^/s ^ 21.7-23.8 GeV. The left (right) plots 
are for theoretical scales fj, = p^/S {p-r/2). Two sets of PDFs (MRST04 and CTEQ6.1M) and three FFs (AKK05, AKK08, 
DSS, from top to bottom) are used. The dashed lines are just indicative for variations of ±30% from the reference at R = 1. 



At large values of p^, the use of the fixed-order perturbation theory is fully justified, since the perturbative 
series is controlled by a small expansion parameter, dsip^)- However, in the typical kinematic range of fixed-target 
experiments, where x^ = 2p,j./y/s ^ 0.1, the coefficients of the perturbative expansion are enhanced by extra 
powers of logarithmic terms of the form d^ In "(1 — x^) or d^ In ^~ (1 — x^). Resummation to all orders of such 
"threshold" terms - which appear because the initial partons have just enough energy to produce the hi gh-t ransverse 
momentum parton - have been carried out at next-to- leading logarithmic (NLL) accuracy in [30|, |3l|. These 
studies confirm that accounting for these terms results in a large (approximately p^-independent) enhancement 
of the perturbative cross section for pion production in the range of fixed-target energies of relevance here 
(\/5 ^ 20 GeV). These studies also find that the scale dependence is also reduced at NLL compared to NLO. The 



presently used fixed-order calculations (INCNLO) do not include threshold resummations but their effect in the final 
spectrum is accounted for, in an effective way, by our choice of relatively small theoretical scales, /x/p^, = 1/2-1/3, 
which results in a cross section increase of a factor of ^--^2-3 as compared to the n/p^ = 1/2-2 range used e.g. in [30l[3ll|. 

The two non-pcrturbative inputs of Eq. ([T]) are the parton densities and the fragmentation functions. The former 
are obtained mainly from global-fit analyses of deep-inelastic electron-proton data, the latter from hadron production 
results in e+e^ collisions. The PDFs are known to within ~20% uncertainty 132] in the kinematic range of interest 
here: x^ = VtIv^'^^ ~ 0.2-0.5 at midrapidity. We use here two of the latest standard PDFs available: MRST04 [sj] 
and CTEQ6.1M [34I. For the quark and gluon fragmentation functions into pions, we use and compare three 
parametrizations: the commonly used AKK05 ^34J plus two more recent sets: DSS [33 and AKK08 [391 ■ The 
dominant fragmentation contribution to Eq. ([T]) comes from the large-z domain: (z) = (phadron/Pparton) ~ 0.8 for 
Py > 3 GeV/c at %fs = 22.4 GeV, where the e+e~ fragmentation data used to obtain the FFs are scarce. In 
addition, the gluon-to-pion FF is not well determined by e^e^ annihilation data, as it appears there only at NLO, 
and we explore small fragmentation scales (in particular, when using [ig — VtI"^^ '^^^ away from the kinematical 
regions where the e+e~ fits are performed. All these issues, which were a concern for the older FF sets like 
KKP [37], Kre [3g| or AKK05, have been partially solved with the most recent fits (2^ which include for the first 
time also hadronic data (and error analyses, such as for HKNS [33|) in their global analyses. These new fits cover 
a larger z range and are more sensitive to the gluon fragmentation. As a result, the normalization of the gluon 
fragmentation function into pions is increased by e.g. up to 50% in AKK08 36] with respect to AKK05 [3J| at the Z'^ 
mass scale. This has an obvious impact in the absolute normalization of the predicted pion spectra as we discuss below. 

In Figs. [2] (spectra) and [3] (ratio data/pQCD), the measured pion p-y single inclusive distributions at various 
energies are compared to the corresponding NLO predictions for varying theoretical scales (^ = Py/3 andpy/2), PDFs 
(MRST04 and CTEQ6.1M) and FFs (AKK05, AKK08 and DSS). In general, the calculations tend to underpredict 
the measured cross sections. The overall agreement, in the p^ dependence and absolute normalization, improves 
going from the left (scales /i = VtI"^^ to the right (scales [i = Pt/'^) and when using MRST instead of CTEQ. The 
MRST04 parametrization results in a cross section 25% larger than using CTEQ6.1M in the range^ p^ = 3-6 GeV/c. 
Such a difference in the resulting cross sections is larger than expected from error analysis within a single PDF set. 
The AKK08 and DSS fragmentation functions reproduce better the data than the AKK05 ones. The overall trend is 
consistent with MRST04 and AKK08/DSS predicting a higher pion yield than CTEQ6.1 and AKK05 in the kinematic 
range of interest here. In any case, the data-theory agreement at fixed-target energies for high-p^ pions is clearly 
better than for prompt-photon production, where the measured E706 yield at -y/s = 31.6-38.6 GeV appears to be two 
to three times larger than the corresponding INCNLO predictions |40l |. 

IV. INCLUSIVE PION SPECTRA IN p-p COLLISIONS AT ^ = 22.4 GEV: A PRACTICAL 

PARAMETRIZATION 

After verifying that the fixed-order pQCD calculations can reproduce relatively well the existing high-p^ pion data 
at fixed-target energies, the second motivation of this study is to provide a practical parametrization of the p-p pion 
spectrum at ^/s = 22.4 GeV to be used as reference baseline for high p^ tt" production in A-A collisions at the 
same cm. energy, where no proton-proton data has been yet measured at RHIC [ijj. We discuss here the method 
followed to obtain a fit from the existing experimental data sets after rescaling them to a common center-of-mass 
energy making use of the NLO predictions. 

A. Center-of-mass energy rescaling 

The existing data sets (Table H]) cover the range of cm. energies from 21.7 GeV to 23.8 GeV. Although at low-p^ 
(below ^ 2 GeV/c), the small differences in yfs result into negligible variations of the soft pion yield and all spectra 
agree well (see Fig. [T]), at high p^ - as one approaches the kinematical limit - a couple of GeV of extra cm. energy 
available can result into a significant change in the parton-parton cross sections. For instance, as can be seen in Fig.|4l 



^ However, closer to the kinematical limit, above 8 GeV/c, the trend changes rapidly and the CTEQ6.1M fit overshoots the MRST04 one 
by up to 40%, indicating the large current uncertainty of the gluon and sea-quark densities at high values of x. 



at p^ =5 GeV/c, going from ^/s = 22 A GeV up (down) to 23.8 GeV (21.8 GeV) results in an increase (decrease) of 
the cross section by a factor of ^ 60% {^ 30%). 
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FIG. 4: Differential vr" cross sections in p-p collisions predicted by NLO pQCD calculations with scales ^ — p^/2, CTEQ6.1M 
parton-distribution-functions, and AKK08 fragmentation functions at ^/s = 21.7, 22.4 and 23.8 GeV. 

Although, as seen in the previous section, there are relatively large uncertainties in the NLO predictions for the 
absolute cross sections, most of these uncertainties cancel out when taking ratios of the predicted perturbative yields 
at different, yet close, cm. energies. In particular, the (large) scale dependence is completely removed. One can, 
thus, rescale all experimental data points measured at a given -/s = X GeV (in the range 21.8 - 23.8 GeV) to a 
common -y/s = 22.4 GeV value via the following prescription: 

rfacxp (Vi = 22.4 GeV) ^ Ua^^^/dp^ (yi ^ 22.4 GeV) \ ^ daexp {^s ^ X GeV) 
dp^ \ da^^^/dp^ (y^ = X GeV) / dp^ 

The pQCD cross-sections are computed in the range p^ « 3-10 GeV/c for the 4 energies under consideration and the 
ratio over the predictions at y/s = 22.4 GeV is fitted to a polynomial form of order 2 or 4. Obviously, to minimize 
the theoretical uncertainties, both the denominator and numerator of the NLO "rescaling factor" (the expression in 
parentheses in Eq. ^) need to be computed using consistently the same PDFs, FFs and scales. The scaling factors 
provided here are obtained averaging over various different choices of these ingredients. The resulting scaling factors 
differ, in any case, as expected by a very small factor ±5%, well covered within the experimental uncertainties alone. 
The functional form of the rescaling factor is chosen so that the correction is zero at p^ =0 GeV/c, so as to obtain a 
smooth extrapolation in the low-p^ region. In any case, below p^ « 1 GeV/c, the correction is (well) below ^5% and, 
so, the experimental low-p^ points are virtually unmodified as they should be by applying this rescaling procedure. 
The final correction functions are shown as a function of p^ in Fig. [5l 

In order to better estimate the uncertainty of the rescaling factors computed theoretically, the energy rescaling 
has also been determined a posteriori^ assuming that the invariant production cross section is a scaling function of 
x^ = 2p^/^/s: 

Ei^^^^i.^,).{±.)'n.^i (3) 

as it should be in perturbative QCD. Taking for the function F the final parametrization discussed in the next Section, 
F[x^) = (22.4 GeYYf{x^ x [11.2 GeV]), the rescaling factor is computed using Eq. ([3]). The difference between this 
empirical estimate and the theoretical rescaling factors, roughly 10%, is assigned as the uncertainty of the present 
energy rescaling procedure. 
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FIG. 5: Rescaling correction factors of the pion cross-sections at cm. energies 21.7-23.8 GeV to a common y^ 
value, as a function of p^ , obtained from the ratio of NLO calculations given by Eq. ^. 



22.4 GeV 



B. Global fit of the rescaled pion Pj, spectra at ^/s — 22.4 GeV 



By applying the appropriate energy correction factors discussed in the previous section to all the experimental data 
sets, we obtain a new set of data points which approximates better the expected 7r° spectrum at y/s — 22 A GeV. The 



experimental spectrum 



Ed"" 



d'^p 



is fitted to the following empirical 4-parameter functional form: 



I{PTAPi\i=0..i) 



Po 



(4) 



Such a formula interpolates well between the low-p^ exponential shape and the high-p^ power-law while fulfilling the 
requirement of being zero at the kinematical limit {p™^^ — 11.2 GeV/c, fixed in the fit). The po parameter gives the 
cross-section at zero pr^ , pi indicates the transition value from soft to hard production, and the p2 and p^ exponents 
characterize the power-law and end of phase-space ranges. After rejecting two data sets which are not consistent with 
the rest of spectra (see below), we obtain a final set of n^at — 194 data points fitted with Eq. ([4]). The resulting fit is 
shown in Fig. [Sj The parameters are obtained from the minimization of the x^ function 



Hfe}) = E 



j=i 



d'^p 



y=o 



(PT,)- f{pT,,{Pi}) 



(5) 



where aj is the statistical and systematic error of point j added in quadrature. The error of the parameters pi are 
given from a deviation of Ax^ from its minimum: 



XH{P^ + SP^})-XH{P^})^^X^ 



(6) 



Although a usual choice is A^^ = 1, we shall conservatively allow for a larger variation of the fit parameters assuming 
Ax^ — 50 in what follows^, similarly to what is done in global fit analyses of parton densities or fragmentation 
functions (see e.g. (35|, |4l|). From this procedure, the corresponding parameters are: 



Po = 176.3 ± 69.7 [mb GeV 

Pi = 2.38 ±1.19 [GeV/c] 

P2 = -16.13±7.21 

P3 = 6.94 ±5.64 

xVndf = 208.2/190 



-'c'] 



(7) 



This would correspond to an increase of 25% of X^in/^'^^ with ndf ~ 200. 



with an important correlation between parameters and errors, as indicated by the large non-diagonal terms of the 
covariance (error) matrix, Vij: 



V^J^ 



1.000 -0.725 

-0.725 1.000 

-0.603 0.981 

0.394 -0.862 



-0.603 0.394 

0.981 -0.862 

1.000 -0.940 

-0.940 1.000 



Note that at low p^, this fit is consistent with an exponentially decreasing function with inverse slope —pi/p2 = 
148 ± 16 MeV. The scale pi, which naively separates soft from hard dynamics, has a sensible value ~ 2-3 GeV. 
Finally, the negative power slope, p2 — —16.13, is found to be larger in absolute value than p2 — —10 obtained at 
\/s = 200 GeV [441 . This is expected from the steeper dependence of parton densities and fragmentation functions 
probed at higher x and z, respectively, at lower ^/s. The uncertainty A/ of the fit is given by 



A/ 



3 

E 

ij=0 



9/ 

dpt 



V,, 



dPj 



1/2 



(8) 



The relative uncertainty of the parametrization. A///, spans the range from ^^15% at low Pj, < 2 GeV/c up to 25% 
(40%) at p^ —A GeV/c (5 GeV/c) in the range covered by the RHIC measurements [i3|. At higher p^'s the fit is 
completely unconstrained due to the lack of data, and its uncertainty is very large. 
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FIG. 6: Compilation of all pion transverse spectra measured in p-p collisions at y^ = 21.7 - 23.8 GeV, rescaled to a common 
•^s = 22.4 GeV energy, as discussed in the text, and fitted to Eq. Q, with the parameters ^. 



Figure [7] shows the ratio of all data sets compiled and rescaled in this work over the fit Eq. ^ with the parameters 
quoted above. All data sets - but Carey76 [lj| at ^/s = 23 GeV and Eggert75 [20| at ^/s — 23.6 GeV which have a 
shape and absolute normalization inconsistent with the rest of measurements and have not been included in the final 
global analysis - show a rather good agreement with the proposed parametrization, as also indicated by x^/ndf ~ 1.1. 
We note that our empirical fit includes also the low-p^ range, not amenable to perturbative analysis, since we want to 
provide a (potentially useful) p-p reference parametrization in the whole range covered by the nucleus-nucleus data. 
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FIG. 7: Ratio of all pion transverse spectra measured in p-p collisions rescaled to -^s = 22.4 GeV over the fit given by Eq. Q 
with the parameters (0 ■ The yellow band represents the uncertainty assigned to the parametrization, as given by Eq. ^ . The 
[carey76] and [eggert75] points at \/s = 23.6,23. GeV respectively il4i. i2Q] have been excluded from the global fit. 

V. SUMMARY 

We have compared the available high-p^ pion spectra measured in proton-proton collisions in the range \/s = 21.7 

- 23.8 GeV (CERN-ISR collider and CERN and FNAL fixed-target) to next-to-leading order pQCD calculations 
with recent parton distribution functions (PDFs) and fragmentation functions (FFs). A choice of the theoretical 
(factorization, fragmentation and normalization) scales between p^/3 and Pj,/2 reproduces well the magnitude and 
shape of the experimental data. CTEQ6.1 and MRST04 parton densities yield results different by up to 25%. 
Second-generation parton-to-pion fragmentation functions (FFs) with updated constraints on the gluon and large-z 
fragmentation region such as DSS or AKK08, improve the agreement of the data with the calculations compared to 
older FF parametrizations. 

A baseline nucleon-nucleon reference p^ -distribution for inclusive tt" production at ^/s = 22.4 GeV has been 
determined from a global fit analysis of the available data. The measured (high-p^ ) data sets have been rescaled at a 
common cm. energy making use of the predicted NLO pQCD yields at the various ^/s. The resulting parametrization 
is consistent within ±15% and ±40% systematic uncertainty with the rescaled 7r° and n^ measurements at low (p^ < 
2 GeV/c) and moderate {pj. ~ 5 GeV/c) transverse momentum. Such a reference - Eq. (|4]) with fit parameters ([7]) 

- can be used in order to obtain the nuclear modification factor of high-p^ pion production in A-A collisions at 
y'S^ = 22.4 GeV measured at RHIC. 
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